Periodizing quasicrystals: Anomalous diffusion in quasiperiodic systems 
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We introduce a construction to embed a quasiperiodic lattice of obstacles into a single unit cell of a higher- 
dimensional space, with periodic boundary conditions. This construction transparently shows the existence of 
channels in these systems, in which particles may travel without colliding, up to a critical obstacle radius. It 
provides a simple and efficient algorithm for numerical simulation of dynamics in quasiperiodic structures, as 
well as giving a natural notion of uniform distribution (measure) and averages. As an application, we simulate 
diffusion in a two-dimensional quasicrystal, finding three different regimes, in particular atypical weak super- 
diffusion in the presence of channels, and sub-diffusion when obstacles overlap. 



PACS numbers: 61.44.Br, 66.30.je, 05.60.Cd, 05.45.Pq 

Quasicrystals are produced by cooling from a melt at a 
rate intermediate between that of periodic crystals (slow) and 
glasses (fast), and have a degree of order which is interme- 
diate between the two, being neither periodic nor random 
0JEI- The transport properties of quasicrystalline materials 
are of particular interest, with diffusion having been exten- 
sively studied experimentally GHZL being of importance for 
their production and technological applications (§)• This is 
related to other transport properties, such as heat conductivity 
and electronic transport [5]; furthermore, the measured ther- 
mopower in quasicrystals is due to electron diffusion (9) . 

To understand such transport properties, it is useful to study 
simple models, in particular the Lorentz gas (LG) fTOlL which 
consists of an array of fixed obstacles in R n with which 
moving particles undergo elastic collisions; The geometry 
in which the obstacles are arranged in such billiard models 
strongly influences these transport properties. 

LGs with a periodic geometry have been extensively stud- 
ied (TTHT5) . Usually, normal diffusion is expected, for 
which the mean-squared displacement has asymptotic behav- 
ior (Ax(t) 2 } ~ t when t — » where Ax(t) := x(t) —x(0) is 
the displacement of a particle at time t, and (•) denotes an av- 
erage over uniform initial conditions. However, a key role is 
played by the presence or absence of channels through which 
particles may travel unimpeded: in the presence of "princi- 
pal horizons", i.e., channels of the highest possible dimension 
n—l (15] [TSJ, there is instead weak super-diffusion with a 
logarithmic correction, (Ax(t) 2 } ~ tint fTsHTSl . 

When the obstacles are rather arranged randomly, normal 
diffusion occurs if the obstacles do not overlap (T9) ; when 
overlaps are allowed, sub-diffusion, with (Ax(t) 2 ) ~ t a and an 
exponent a < 1 , is expected, due to the presence of arbitrarily 
large traps near the percolation threshold (201I2T) . 

It is then natural to ask what type of diffusion is expected 
for the intermediate case of a quasiperiodic LG, as was raised 
in ref. [ 22 ] for the special case of obstacles arranged in a Pen- 
rose tiling. However, we are not aware of any analytical or 
numerical results on such systems, except for a rather non- 
physical one-dimensional model (23). Indeed, the numerical 
study of quasiperiodic and random systems is known to be 



challenging, due to the absence of periodic boundary condi- 
tions, which in the periodic case allow us to restrict attention 
to a single unit cell. For non-periodic systems, it is often nec- 
essary to generate at the large structures at the outset, leading 
to inefficient algorithms (23). 

In this Letter, we introduce a construction to study 
quasiperiodic systems by embedding them in a periodic sys- 
tem of higher dimension, which works by reversing the pro- 
jection method used to construct quasiperiodic lattices O 
|25) . This construction solves the above problem for quasiperi- 
odic systems, allowing us to again reduce the system to a sin- 
gle unit cell with periodic boundary conditions, but now in the 
higher-dimensional system. 

The construction has several consequences: Firstly, it 
shows in a simple way the existence of channels in quasiperi- 
odic structures, up to a critical radius of the obstacles. Sec- 
ondly, it gives a natural notion of uniform distribution (mea- 
sure) in the system, and hence of averages, by reducing the in- 
finite system to a finite (compact) one. Finally, it leads directly 
to a simple and efficient algorithm for simulating dynamics in 
quasiperiodic structures. 

By applying this simulation method, we observe three dif- 
ferent diffusive regimes for a two-dimensional quasiperiodic 
LG, including weak super-diffusion in the presence of chan- 
nels, and sub-diffusion when the obstacles are overlapping. 

Projection method:- The projection method O I26U29) 
constructs a quasiperiodic lattice in a subspace E (includ- 
ing the origin) of a Euclidean space R n by projecting a sub- 
set of the vertices of a hypercubic lattice L, consisting of 
points with integer coordinates in R n , onto E\ for example, 
the well-known Penrose tiling can be obtained by projecting a 
5-dimensional cubic lattice into two dimensions |24-26|. We 
denote by m the dimension of E, with m <n, and by E± the 
orthogonal complement of E, of dimension n — m, such that 
E ® E± = R n . In the following, for simplicity, we refer to 
orthogonal projections simply as "projections". 

In one version of the projection method (26) , we choose 
which lattice points to project by considering the Voronoi re- 
gion of a lattice point p of the cubic lattice L, i.e., a cube cen- 
tered at p, and projecting onto E exactly those lattice points 
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FIG. 1: (Color online) Two versions of the projection method: (a) 
Using Voronoi regions, shown as shaded (green) boxes. The arrows 
indicate the lattice points in L C R 2 which are projected onto the 
line E with irrational slope, (b) Canonical version, projecting lattice 
points lying inside a strip. 
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FIG. 2: Embedding of a 1 -dimensional quasiperiodic lattice into a 
two-dimensional periodic billiard, (a) Construction of the set P by 
projection and "thickening". The parallel diagonal lines (red online) 
show different "slices" E y . (b) The billiard obstacle K, after peri- 
odizing P, shown as the oblique (red) bar. Examples of short and 
long paths for the resulting billiard dynamics are shown. 



p whose Voronoi regions intersect E\ see fig.[TJa). This gives 
a set L|| of points in the subspace E, which is a quasiperiodic 
lattice if E is totally irrational, i.e., the angle made with each 
lattice direction is an irrational multiple of n l26li . 

The Voronoi projection method in fact produces the same 
quasilattice as the following canonical projection method 
l26l : consider a cube whose side length is the lattice spac- 
ing, centred at the origin, and project it onto the orthogonal 
subspace E±, giving a set W. The lattice points in L which lie 
inside W x E are projected onto E\ see fig.[TJb). 

Construction of embedding:- The quasiperiodic LG con- 
sists of balls of dimension m in E centred at each point of L\\ . 
We will, in a sense, reverse the projection method, to con- 
struct a billiard obstacle K inside a unit cube C of R n with 
periodic boundary conditions, designed specifically such that 
trajectories of the billiard dynamics parallel to E inside C au- 
tomatically give trajectories of the quasiperiodic LG when the 
dynamics is "unfolded" to the whole of R n - i.e., by following 
the dynamics using periodic boundary conditions in one cell, 
but keeping track of which cell has been reached l30ll . 

We start by considering free motion of particles with arbi- 
trary initial positions inside C, but whose velocities are con- 
strained to move parallel to E. When they cross the periodic 
boundaries, they keep the same velocity, but jump to a differ- 
ent translated copy, or "slice", E y := E + v, parallel to E but 
translated by a vector v which is orthogonal to E, i.e., such 
that v is in E± . The periodic boundary conditions can thus be 
thought of as allowing us to "wrap" the whole of the plane E 
in R n inside one single unit cell C. 

Since E is totally irrational, a trajectory emanating from a 
generic initial condition in the cube C will fill C densely Oft ; 
similarly, the collection of parallel slices E y which are visited 
by a given trajectory also fill C densely. We will construct the 
billiard obstacle K inside C such that the dynamics retains this 
denseness property. 



Since E is wrapped inside C, we can think of C as corre- 
sponding to each of the Voronoi cells in R n which intersect E. 
The projection method then requires that we project the origin 
onto each slice E y which intersects C. We denote by W the 
set of all such points; it turns out to be exactly the projected 
set used in the canonical projection method. Note that W is a 
subset of full dimension inside E±, which can be constructed 
in general by projecting the vertices of the cube C onto E±_ 
and taking the convex hull of the resulting points. 

To complete the construction, we now define a set P as the 
result of projecting not just a single point at the origin, but 
rather the whole billiard obstacle in E, an m-dimensional ball 
B of radius r centred at the origin, so that P := W x B. How- 
ever, some parts of P constructed in this way fall outside C; see 
fig.[2ja). We finish by periodizing P via the periodic boundary 
conditions, giving a new set K, also of dimension n, which is 
the billiard obstacle that we finally use; see fig. |2jb). Since 
the billiard particles move parallel to E, and the boundary of 
the obstacle K is cylindrical, with axis perpendicular to E, the 
particles will always move parallel to E in IR W . When the re- 
sulting billiard trajectory is unfolded to R n , it then gives rise 
to a billiard trajectory in the quasiperiodic LG. 

As a simple example, consider a one-dimensional 
quasiperiodic billiard with n = 2 and m = 1, so that C is a 
square, E is a straight line through the origin with irrational 
slope a, and the obstacle B becomes a line segment of length 
r. The set K then consists of three thickened line segments 
with slope —1/oc, as shown in figure [2jb). Since K divides 
the square completely into two parts, there is, as expected, no 
diffusion in this case: any given trajectory remains confined, 
bouncing between two neighboring obstacles. However, this 
point of view easily allows us to calculate, for example, the 
relative probabilities of long and short paths starting from ran- 
dom initial conditions, by calculating the size of the respective 
areas. 
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FIG. 3: QLG en 2D embedded in a 3D cube. The three pieces of 
cylinder form the billiard obstacle K, and the billiard dynamics takes 
place in the planes E y orthogonal to K. There are two planar channels 
(not shown explicitly) lying along the cube's faces, and an additional 
channel shown at an angle of J which does not intersect K for this 
choice of obstacle radius. 



2D quasiperiodic Lorentz gas:- We now turn to the 
simplest non-trivial application of this construction, a 2D 
quasiperiodic LG formed by projecting a 3D simple cubic lat- 
tice with cubic unit cell C. The subspace E is a totally irra- 
tional plane through the origin, and the obstacle B is a disc of 
radius r. The resulting 3D periodic billiard model is shown 
in figure [3| the corresponding billiard obstacle K consists of 
three segments of a cylinder, in analogy to the previous sec- 
tion. The two parts of the cylinder which come from periodiz- 
ing W are capped by those planes E y which pass through the 
respective vertices of C. 

The construction provides a natural way to interpret a uni- 
form distribution of initial conditions in the phase space of 
the quasiperiodic LG: positions are uniform in the cube C, 
and velocities have unit speed and point in uniform directions 
parallel to E. We conjecture that the dynamics is ergodic, i.e., 
from almost any starting point outside the billiard obstacle K 
(except for a set of measure zero), the corresponding trajec- 
tory fills densely and uniformly this phase space. Averages 
are then taken with respect to this uniform distribution. 

Channels:- Motivated by the results on periodic Lorentz 
gases, we ask if it is possible that there exist channels through 
the quasiperiodic LG structure, which could lead to super- 
diffusive behavior |[T5l[T6l . The absence of periodicity leads 
to the expectation that this is not possible, as suggested 
in ref. l22l . However, our construction in fact shows that 
quasiperiodic LGs do have channels when the obstacle radius 
is sufficiently small. 



Indeed, suppose that the radius r of the obstacle S, and 
hence of the cylinder K, is small enough that none of the 
three cylindrical pieces (figure[3} intersect the remaining faces 
of the cube (i.e., those not intersected by their axes). Then 
in the 3D billiard, there are "planar channels" n [15 ] lying 
along each of these faces. Restricting the dynamics to be 
parallel to the plane E, we see that each planar channel n 
in the 3D periodic billiard induces a rectangular channel in 
the 2D quasiperiodic LG, given by the intersection of E with 
n. Elsewhere we will show that this result generalizes to 
quasiperiodic LGs in m > 2 dimensions, which have (m — 1)- 
dimensional "principal horizons" fl6l . Note that as the radius 
decreases, additional channels may appear, as is also the case 
in the periodic LG (Tt} [T8]|. For example, fig. [3] shows such 
additional planar channels oriented at an angle 71 / 4. 

Numerical results:- The above construction provides a 
direct and simple simulation method for quasiperiodic sys- 
tems, by simulating directly billiard dynamics in the higher- 
dimensional system with obstacle K and periodic boundary 
conditions. We applied this to the 2D quasiperiodic LG to in- 
vestigate its diffusive properties as a function of the radius r 
of the billiard obstacles. 

In our simulations, we place 10 6 particles distributed with 
uniform positions in the three-dimensional cube C, and veloc- 
ities with unit speed distributed with uniform directions par- 
allel to the subspace E, and calculate the mean-squared dis- 
placement (MSD) (x(t) 2 ) as a function of time. The subspace 
E± was taken aligned along the unit vector ( ? 0+2 > ^0+2 )' 
where (j) := 1+ 2 ^ is the golden ratio. 

Figure [4] shows (x(t) 2 )/t, in order to emphasize deviations 
from normal diffusion. We find that for small radii, this quan- 
tity increases in time, so that super-diffusion is present. How- 
ever, the correction does not appear to be a simple logarithmic 
factor, as for periodic LGs with channels, but rather seems to 
be fit well by expressions of the form (x(t) 2 } ~ t(lnt) s , with 
an exponent 8 < 1 which depends on the radius r. 

This super-diffusive behavior occurs due to the presence of 
channels. As r increases, the channels are blocked one by 
one, as the billiard obstacles expand sufficiently to intersect 
the planes defining the channels. Note that once the billiard 
obstacles are large enough to cross the faces of the cube, their 
periodic images must also be taken into account as additional 
obstacles in the simulations; for this reason, the algorithm is 
most efficient when the obstacles are small enough for this not 
to occur, which is exactly the most difficult case to treat using 
standard methods, for which very large quasiperiodic lattices 
must be constructed. 

At a certain critical radius around r ~ 0.21 (discussed in 
more detail elsewhere), all of the planar channels become 
blocked, but the billiard obstacles still do not touch each other. 
This is the equivalent of the "finite horizon" regime in periodic 
LGs lfl31[T6l . and the simulation results for the quasiperiodic 
LG indeed appear to show normal diffusion, with (x(t) 2 ~ t, 
in this case. However, as r increases further, there is another 
critical radius around r ~ 0.30, at which the billiard obstacles 
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FIG. 4: (x(t) 2 )/t as a function of t for different radii r, shown in (a) 
semi-logarithmic and (b) logarithmic scale, with a zoom in the inset. 




FIG. 5: Representative trajectories in the 2D quasiperiodic LG; the 
gray scale (color online) indicates time, as shown in the respective 
color bars. Blank, circular regions correspond to the billiard obsta- 
cles; the positions of the obstacles are identical in both figures, (a) 
r = 0.2; several initial conditions are shown, (b) r = 0.425; single 
initial condition. 

touch and then begin to overlap one another. In this regime 
we find sw/?-diffusion (highlighted in fig. |4jb)) as in the ran- 
dom, overlapping LG, until the final critical radius r ~ 0.43 
is reached, when all particles become localized in bounded 
regions of space and there is no longer diffusion. 

Figure [5ja) shows several representative trajectories for 



small r in the super-diffusive regime, highlighting the chan- 
nels within the quasiperiodic LG which give rise to the super- 
diffusive behavior. For the radius r shown, there are three 
directions in which channels occur, corresponding to those in 
fig. [3]- two along two perpendicular faces of the cube, and 
one at an angle of |. Figure J5^b) shows a single trajectory for 
large r, close to the percolation threshold at which diffusion 
ceases. Narrow bottlenecks between cavities of different sizes 
are visible, which explains the observed sub-diffusive behav- 
ior, in analogy to the overlapping random Lorentz gas. 

In summary, we have introduced a construction to study 
quasiperiodic lattices by embedding them into a periodic 
unit cell in a higher dimensional space. This point of view 
transparently shows that there quasiperiodic Lorentz gases 
have channels for small obstacle radius, and provides a sim- 
ple, direct and efficient simulation method for dynamics in 
quasiperiodic structures. Numerical simulations exhibit three 
diffusive regimes, including super- and sub-diffusion. We ex- 
pect that this construction can be profitably applied to simulate 
more complicated quasiperiodic systems. 
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